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Abstract 

We studied 101 flux emergence events ranging from small ephemeral regions to large emerging flux 
regions which were observed with Hinode Solar Optical Telescope filtergram. We investigated how the 
total magnetic flux of the emergence event controls the nature of emergence. To determine the modes of 
emergences, horizontal velocity fields of global motion of the magnetic patches in the flux emerging sites 
were measured by the local correlation tracking. Between two main polarities of the large emerging flux 
regions with more than around 2 x 10 19 Mx, there were the converging flows of anti-polarity magnetic 
patches. On the other hand, small ephemeral regions showed no converging flow but simple diverging 
pattern. When we looked into the detailed features in the emerging sites, irrespective of the total flux 
and the spatial size, all the emergence events were observed to consist of single or multiple elementary 
emergence unit(s). The typical size of unitary emergence is 4 Mm and consistent with the simulation 
results. 

From the statistical study of the flux emergence events, the maximum spatial distance between two 
main polarities, the magnetic flux growth rate and the mean separation speed were found to follow the 
power-law functions of the total magnetic flux with the indices of 0.27, 0.57, and -0.16, respectively. From 
the discussion on the observed power-law relations, we got a physical view of solar flux emergence that 
emerging magnetic fields float and evolve balancing to the surrounding turbulent atmosphere. 

Key words: Sun: magnetic fields — Sun: emerging flux — Sun: photosphere — Sun: chromosphere 
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1. Introduction 

The sites where the subsurface magnetic flux tubes 
emerge on solar surface are called emerging flux regions 
(Bruzek 1969; Zirin 1972). The typical simple emerging 
flux region has one main pair of the opposite magnetic 
concentrations at the both ends of the emerging site. The 
opposite magnetic concentrations of the main pair move 
away from each other with the speed of 1-2 km -1 (Zwaan 
et al. 1985; Brants 1985) in the developing phase of the re- 
gion. Inside the developing emerging site there are many 
magnetic flux tube emerging on the photosphere. The 
newly emerged magnetic flux tubes at the emerging site 
appear as the dark granular lanes (Loughhead and Bray 
1961; Brants and Steenbeek 1985) on the photosphere. At 
the both ends of the dark granular lane, there are the mag- 
netic concentrations which are called the footpoints and 
observed as the bright points in G band image (Otsuji et 
al. 2007). With Ha the emerged flux tubes are observed as 
dark arch filaments (Bruzek 1967). The lifetime of an arch 
filament is 10-30 minutes (Bruzek 1967; Chou and Zirin 
1988). The rise velocity of arch filaments is 10-15 km s _1 
(Bruzek 1969; Chou and Zirin 1988). Otsuji et al. (2010) 
found the deceleration of the apex of the small-scale arch 
filaments in the chromosphere. 

Emerging flux regions show a variety of size, lifetime, to- 
tal magnetic flux and field strength. Especially the small 
emerging flux regions are called ephemeral active region 
(EAR; Harvey and Martin 1973). For convenience, in the 



following page we define larger (i.e. non-EAR) emerging 
flux regions as EFRs. EFRs are produced by fairly large- 
scale flux emergence. They have a pair or more complex 
group of sunspots with definite penumbrae. The typical 
size of EFRs is more than 30 Mm (Bruzek 1967). EFRs 
show their emergence activities for several days and exist 
on the solar surface for a few month at the maximum. The 
total flux in an EFR increases with the rate of 10 20 Mx 
hr _1 and reaches to the order of 10 20 -10 22 Mx (Zwaan 
1987). In the fairly developed EFR the field strength 
of main spot is around 3,000 Gauss (Brants and Zwaan 
1982). 

EARs have simple bipolar configuration. They have no 
penumbra in the sunspots. The typical size varies from 
5 Mm to 30 Mm (Harvey and Martin 1973; Harvey et 
al. 1975; Hagennar 2001; Otsuji et al. 2007). EARs have 
short life time of hours or one day. The total flux in an 
EAR is up to 10 20 Mx with the increase rate of 10 19 Mx 
hr . The magnetic field strength of main spot is from a 
few times 100 Gauss (Martin 1988) to 2,000 Gauss (Brants 
and Zwaan 1982). 

Various simulation studies on flux emergence have 
been performed by many researchers (Shibata et al. 
1989; Matsumoto and Shibata 1992; Fan 2001). They 
showed the simple bipolar emergence simulation which 
corresponds to the observation result. Matsumoto et al. 
(1993) and Magara and Longcopc (2001) simulated the 
three-dimensional magnctohydrodynamics (MHD) of the 
emerging magnetic flux. Nozawa (2005) and Murray et al. 
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(2006) also performed the MHD simulations of flux emer- 
gence with the sheared or twisted flux tube. They found 
that the flux tube with shear or twist emerges faster than 
that without any shear and twist. 

Nozawa et al. (1992) performed the MHD simulation of 
flux emergence in a sheet geometry. The initial stable flux 
sheet in the convective zone was perturbed with various 
wavelengths, which correspond to the convective motion. 
They found that irrespective of the wavelength of initial 
perturbation, a finite "most unstable wavelength" is ex- 
cited. This wavelength (2-4 Mm) is inherent in the Parker 
instability (Parker 1966). As a result, the flux sheet is un- 
dulated and the apexes of the convex field line (il-loops) 
appear consecutively on the photosphere. Some of dipped 
field lines (U-loops) also emerge to form the regions called 
"bald patches (Titov et al. 1993)" . 

On the other hand, observational study on a large EFR 
and bald patches was performed by Pariat ct al. (2004). 
They observed a fairly large (<~ 30 Mm) EFR with the 
magnetogram and found that the emerged field lines un- 
dulate vertically. They revealed that there are many bald 
patches between the main spots. These results confirmed 
that the emerging flux tube does not rise altogether at a 
time, but each il-loop component rises individually. They 
proposed this model as "resistive emergence model" . The 
distance between two consecutive bald patches is in the 
range of 2-6 Mm, which is consistent with a theoretical 
argument on the flattening of emerging magnetic field just 
below the surface and its critical emergence length, first 
presented by Magara (2001). 

Recently Isobe et al. (2007) further developed the sim- 
ulation performed by Nozawa et al. (1992) and obtained 
the result in which the undulated field line caused re- 
connections with neighboring 17-loops and created larger 
loops. Archontis and Hood (2009) performed the three- 
dimensional MHD simulation of the emergence of undulat- 
ing fieldlines. These reconnection events are interpreted as 
a sources of Ellerman bombs (Ellerman 1917; Kurokawa et 
al. 1982; Kitai 1983; Matsumoto ct al. 2008a; Matsumoto 
ct al. 2008b; Watanabe et al. 2008). 

The resistive emergence model is applicable to the 
large EFR. Furthermore, recently Otsuji et al. (2007) 
found bald patches inside a small-scale EAR (~ 5 Mm) 
using Solar Optical Telescope (SOT; Ichimoto ct al. 
2004; Tsuneta et al. 2007; Suematsu et al. 2007; Shimizu 
et al. 2007) aboard Hinode (Kosugi et al. 2007). However, 
Centeno et al. (2007) and Gugliclmino et al. (2008) showed 
small EARs (2 Mm and 6 Mm, respectively) without un- 
dulated magnetic field. 

As stated above, the criteria of the bald patches for- 
mation are still uncertain. Furthermore, in the latest 
simulation, the footpoints of emerged flux loops showed 
converging motion toward the bald patches on the pho- 
tosphere (Cheung et al. 2010). Although this converging 
motion was observed in preceding studies (Strous ct al. 
1996; Strous and Zwaan 1999; Bernasconi 2002; Cheung 
et al. 2008), statistical analysis on that phenomenon with 
respect to the size and magnetic characteristics of the flux 
emergences has not been done yet. 



To clarify the criteria of forming bald patch and con- 
verging flow, wc performed a statistical study about the 
nature of magnetic flux emergences using SOT. The flux 
emergence phenomena from small EARs to large EFRs ob- 
served by SOT were investigated on their morphological 
and magnetical characteristics. Furthermore, we derived 
the the relations between the total magnetic flux and the 
maximum spatial size, the flux growth rate and the mean 
separation speed of the emergence event to clarify how 
the total flux amount controls the entire evolution of the 
emergence. 

2. Observation and Data Reduction 

2.1. Observation and Data Selection 

Hinode satellite has observed solar surface for over 4 
years with SOT. SOT has Broadband Filter Imager (BFI) 
and Narrowband Filter Imager (NFI). Call H (3968.5 A) 
filtcrgrams were taken by BFI with the bandpass of 3 A. 
Fel (6302 A) and Nal D (5896 A) polarimetric data in 
solar photosphere and chromosphere were observed with 
NFI. To search for the emerging flux phenomena, we used 
Hinode daily quicklook movies 1 . Among the possible can- 
didates, we selected 101 emerging flux phenomena accord- 
ing to the criteria as follows: (1) conspicuous presence of 
separating bright points in Call H image and/or oppo- 
site polarities in Fel or Nal D image, (2) location fairly 
inside the solar limb (a = arcsin^/i?) < 70°, where r is 
the distance from solar disk center to the location of the 
event and R is the solar radius) (3) observational time 
span longer than 1 hour, and (4) observational cadence 
higher than 10 minutes. The period in which we studied 
is between 2006 November 26 to 2010 August 23. The 
detailed observational information of the emerging flux is 
available in the associated electronic tables 2 . 

2.2. Data Reduction 

In this section we give the description of the data reduc- 
tion with the data of EFR 20061201 (# 2 in the electronic 
table) as an example. First, dark-current subtraction and 
flat fielding were performed on the obtained SOT data in 
the standard manner. Then we processed every observed 
data as described below. 

First, we used 50770-Michclson Doppler Imager (MDI; 
Scherrer et al. 1995) magnetogram data to calibrate SOT 
polarimetric data. We compensated the differential ro- 
tation of two consecutive MDI data observed before and 
after the SOT observation and interpolated them by time 
to estimate the distribution of magnetic field at the time of 
SOT observation. Next we deteriorated SOT polarimet- 
ric data with the spatial resolution of MDI magnetogram 
(2"). Using the deteriorated SOT polarimetric data and 
MDI magnetogram, we made a scatter plot of SOT polari- 
metric signal (Stokes V/I) to MDI field strength (Figure 
2). The correlation coefficient of SOT V/I to MDI field 

1 http:/ /hinode. nao.ac.jp/QLmovics/ 

Before the publication of this paper, the table is shown in 
http:/ /www. kwasan.kyoto-u.ac.jp/~otsuji/electronic_table.pdf. 
After publication the table will be uploaded in PASJ site. 
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Fig. 1. (a) Pointing-corrected SOT polarimctric image with 
the white box indicating the area of the "top-view image" 
in Figure 4. (b) Deteriorated SOT polarimctric image, (c) 
Reference SOHO MDI image estimated by linear interpola- 
tion. 

strength was 0.95. We performed linear fitting on the scat- 
ter plot and obtained the conversion equation from SOT 
V/I to the photospheric field strength B p , expressed as 

B p ^B 1 xV/I + B a . (f) 

In the case of Figure 2, the offset value Bq and the scaling 
factor Bi were 4.029 Gauss and 12528 Gauss, respectively. 
Note that the scatter of data points in the figure is mainly 
due to the Doppler effect arising from the satellite's orbital 
motion. With this method we converted the SOT polari- 
metric signals to the magnetic flux densities for all the 
samples. 

In measuring the actual size of the solar features, we 
compensated the projection effect. Then we applied the 
subsonic filter of 3 minutes to the SOT Call H image 
sequence in order to suppress chromospheric oscillatory 
motions (Figure 3b). 



SOT V/I vs. MDI magnetogramw 




SOT V/I 

Fig. 2. Scatter plot of the polarimetric signal V/I of SOT 
and the field strength of MDI magnetogram. 

3. Analysis 

As is stated in the introduction, mutually approach- 
ing anti-polarity patches were observed in undulating re- 
sistive emerging phenomena (Strous and Zwaan 1999). 
We studied the formation process of the converging mo- 
tion between the opposite polarities in EFR and EAR 
by analysing the morphological, dynamical and magnetic 
evolutions of our SOT samples. In this section we in- 
troduce our analysis methods using the data of EFR 
20061201. 

3.1. Morphological Evolution 

The morphological evolution of the magnetic flux emer- 
gence was analysed with two methods. One is the 
method of tracking magnetic elements using local correla- 
tion tracking (LCT), and the other using time-sliced dia- 
gram. 

3.1.1. Local Correlation Tracking 

The local correlation tracking (LCT) method is com- 
monly used to derive the horizontal velocity field 
(November and Simon 1988; Berger et al. 1998; 
Matsumoto and Kitai 2010). For LCT, we used 
f lowmap . pro in SSW of IDL. Flowmap . pro calculates the 
two dimensional vector flowfield by following the subtiles 
in the time series of two dimensional images. 

To examine the motion of the footpoints of flux tubes 
for all over the emerging site, we performed LCT on SOT 
magnetogram data and obtained velocity field of moving 
magnetic elements (Figure 3c). The size of tracking sub- 
tile for LCT was 0"5. To reduce the velocity noise due to 
LCT error, the velocity fields were averaged both spatially 
and temporally over the zones of 1" x 1" and 10 minutes, 
respectively. The standard deviation of the velocity field 
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inside the data cube of 1" x 1" x 10 minutes was ~ 0.1 km 
s- 1 . 

Then we derived divergence of horizontal velocity field 
(Figure 3d). To emphasize global and sustained flows, the 
divergence maps were averaged both spatially and tempo- 
rally. The spatial average was performed with the width 
of 10" for large EFR and 3" for small EAR. The thresh- 
old between large EFRs and small EARs was fixed to be 
40", which is the upper limit size of ephemeral active re- 
gions indicated by Harvey and Martin (1973). For EFR 
20061201, the spatial averaging box size was 10". The 
temporally averaging period was taken as 10 minute for 
all the events. 

3.1.2. Time-sliced Diagram 

To clarify the dynamics of footpoints more quanti- 
tatively, we made time-sliced diagram of SOT magne- 
tograms (Figure 3c). The spatial slit was located parallel 
to the axis of the EFR (shown in Figure 3a). From the 
time-sliced diagram, we derived the maximum distance 
d max between the main spots, the mean separating speed 
(v) of the main spots. Detailed footpoints motions, such 
as mutually approaching anti-polarities, were examined in 
this diagram. 

3.2. Temporal Evolution of Magnetic Field 

To investigate the temporal evolution of emerging mag- 
netic flux, we measured the total flux within the emerging 
site. Positive and negative fluxes were summed up sepa- 
rately. The total flux of the emerging region $ was derived 
by subtracting the fluxes at the initial time. We plotted 
the variation of positive and negative fluxes, respectively 
(Figure 3f). From this plot, we derived the maximum 
amount of unsigned total magnetic flux $ max and the un- 
signed flux growth rate (d$/dt) = $ max /T. The growth 
rate is defined as the total magnetic flux $ max divided by 
the continuously emerging period T. Note that the unbal- 
ance between positive and negative fluxes in the sample 
plot was due to the flow of the positive flux out of the field 
of view. 

4. Results 

First, we introduce the sample results for large EFR 
and small EAR at section 4.1 and 4.2, respectively. Then 
the statistical results are shown in section 4.3. 

4.1. Large EFR 20061201 

In Figure 3a, there are two main spots aligned east- 
west direction. The size of the main spots was 10"-20". 
The field strength of the spots was ±1.8 x 10 3 Gauss at 
the maximum. Although the following spot is missing 
in the Call H image because of the field of view limita- 
tion (Figure 3b), there are sunspots locating at the same 
position with the preceding negative spots in the mag- 
netogram. The velocity field shows prominent outward 
motions of the main spots (Figure 3c). There are pos- 
itive divergence areas inside the main spots (Figure 3d), 
which indicates the flux emergences. From the time-sliced 
diagram (Figure 3e), we can estimate the speed of main 



spots to be ~ 0.3 km s _1 for each, thus the mean sepa- 
rating speed (v) is about 0.6 km s" 1 . At the end of the 
observation period, the distance between two main spots 
increased to around 100". We considered this value as the 
maximum distance d max . The emergence started at 16:00 
UT on 1 December and lasted until the observation end 
at 24:00 UT on 2 December. Thus the active emergence 
period was taken as 32 hours. The maximum amount of 
total magnetic flux were 1.3 x 10 21 Mx for positive polar- 
ity and 3.1 x 10 21 Mx for negative polarity (Figure 3f). As 
the following positive spot flowed out of the field of view, 
measured positive flux was less than that of negative one. 
Thus we took the maximum negative flux as <J> max . Mean 
flux growth rate (d$/dt) for this event was 9.7 x 10 19 Mx 
hr- 1 . 

Let us look at the central part of the region where there 
are many small positive or negative magnetic patches 
(Figure 3a). These patches correspond to the Call H 
bright points (Figure 3b). Although the magnetic patches 
seem to move with apparently random velocities in Figure 
3c, these patches are located in the converging region of 
divergence map (Figure 4d). Thus these magnetic patches 
accumulated and stagnated to the localized area. In fact, 
the time-sliced diagram shows that these magnetic patches 
actually approach to each other with speed of ~ 1 km 
s . So we identify these mutually approaching area as 
a U-loop formation in the EFR. In the following, we de- 
note the area as the "stagnation zone (SZ)", where the 
anti-polarities mutually approach and accumulate. 

4.2. Small EAR 20070213 

We selected a small-scale magnetic emergence event 
which emerged on 13 February 2007 (# 9 in the elec- 
tronic table) as the sample case of EAR. Figure 4a shows 
the magnctogram of the region. Two magnetic concentra- 
tions in the magnetogram correspond to the Ca II H bright 
points in Figure 4b. The field strength of the two mag- 
netic concentrations was ±300 Gauss at the maximum. 
The velocity field of the small EAR shows separative and 
anti-clockwise rotational motion of two magnetic concen- 
trations (Figure 4c). Figure 4d is the divergence map de- 
rived from the velocity field and averaged spatially with 
3". The divergence map shows the positive area at the 
central region of the EAR, which indicates that there is 
no converging flow inside the emerging site. In the time- 
sliced diagram we can see the simple separating motion 
of two main magnetic concentrations and no stagnation 
zone (Figure 4d). The maximum distance between the 
main concentration d max was 10". The mean speed of 
separating motion for main concentrations (v) is about 
1.5 km s _1 . The temporal evolution of total flux shows 
that the maximum amounts of total fluxes are 5.4 x 10 18 
Mx for positive polarity and 4.2 x 10 18 Mx for negative 
polarity (Figure 4f). We took the maximum positive flux 
as $ max for this event. The time-sliced diagram and the 
total fluxes evolution plot indicate that the flux emergence 
began at 03:15 UT and last until the observation end at 
04:30 UT. Thus the active emergence period was 75 min- 
utes and the mean flux growth rate (d$/di) was 4.3 x 10 18 
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Fig. 3. EFR 20061201. (a) SOT Fel magnetogram. (b) SOT Call filtcrgram. (c) Horizontal velocity field of magnetogram derived 
by LCT, averaged with 1" and 1 minutes. Red and blue arrows indicate the velocities of positive and negative magnetic components. 
The threshold field strength of drawing arrows is ±50 Gauss, (d) The divergence map derived from the horizontal velocity field. 
Spatially averaging width is 10". Temporally averaging period is 10 minutes, (e) Time-sliced diagram of SOT magnetogram. The 
both edge of the slit (A and B) in panel (a) correspond to 0" and 120" in the plot, respectively. The vertical white gaps represent 
the observation breaks. The mutually approaching anti-polarities are indicated by black boxes, (f) The evolution of magnetic flux 
of the EFR. 
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Fig. 4. EAR 20070213. (a) SOT Fel magnetogram. (b) SOT Call filtergram. (c) Horizontal velocity field of magnetogram derived 
by LCT, averaged with 1" and 1 minutes. Red and blue arrows indicate the velocities of positive and negative magnetic components. 
The threshold field strength of drawing arrows is ±50 Gauss, (d) The divergence map derived from the horizontal velocity field. 
Spatially averaging width is 3". (e) Time-sliced diagram of SOT magnetogram. The both edge of the slit (A and B) in panel (a) 
correspond to 0" and 20" in the plot, respectively, (f) The evolution of magnetic flux of the EAR. 



Mx hr^ 1 . 

Jf.,3. Statistical Result 

Table 1 in Appendix 1 shows the measured quantities 
of all the samples. If there was no magnetic observation, 
Call H data were used to derive d max and (v), and to 
judge the existence of SZs. There is no data of the total 
flux $ max and the flux growth rate {d$/dt) for the ob- 
servations without the magnetogram. Figure 5 shows the 
statistical characteristics of measured quantities. 
4- 3.1. Existence of SZ 

In Figure 5a, the histogram of the maximum separa- 
tion distance d max indicates that SZs are rarely found in 
small d max regions. The flux emergence phenomena with 
the separation size of more than 25 Mm always have SZs. 
Below the threshold of 25 Mm, there are both EFRs and 
EARs with/without SZs. To clarify the conditions to have 



SZs, we categorized all the events to three groups as follow. 
Group I: maximum separation distance d max > 25 Mm 
and with SZs. Group II: maximum separation distance 
dmax < 25 Mm and with SZs. Group III: maximum sepa- 
ration distance d ma x < 25 Mm and without SZ. Group I, 
II and III are indicated in the Figure 5 by white, gray and 
black bars/squares, respectively. Figure 5b shows the his- 
togram of the maximum fluxes f mlXl where we present the 
maximum flux separately for these three groups. Group I 
has the maximum flux of $ max ~ 10 21 Mx, while group II 
and III have $ max - 10 20 Mx and - 10 19 Mx, respectively. 
Figure 5c shows the relation between the maximum flux 
"JVax and the maximum separation distance d max . In the 
plot the three group I, II and III are clearly separated. 
The maximum separation distance d m ax depends on the 
maximum flux amount $ mal . The scatter plot implies a 
power-law relation of 
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Fig. 5. Statistical characteristics of magnetic flux emergences. Data points are plotted with different symbols according to the 
association with/without the stagnation zone (SZ). (a) Histogram of d max . (b) Histogram of $ max . (c) Scatter plot of <t> ma x and 
<^max- The solid line in the plot indicates the distribution relation of c£ m ax — 7.9 x lO -5 *^. (d) Scatter plot of draax and ^ un it ■ 
(e) Scatter plot of "I> m ax and (d<J>/di). The solid line in the plot indicates the distribution relation of (d^/dt) = 9.6 X 10 7 <t>[} 1 ! 7 { . (f) 
Scatter plot of "3> ma x and (v). The solid line in the plot indicates the distribution relation of (v) = 2.1 X 10 3 $max 16 - The dashed lines 
in the panels represent the deviation (1<t) of the fitting plots. 
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C = 7.9 x 10" 5 $ 



-5^0.27 



(2) 



where d max is in Mm and $ max is in Mx. We can see the 
trend in which the small EARs have low values of maxi- 
mum flux while large EFRs have high values of maximum 
flux. The power-law relation is consistent with the re- 
sult of Hagcnnar (2001), although the index of power was 
0.18 instead of 0.27. From the result of categorization, we 
found that the SZ features are associated with the mag- 
netic flux emergence of more than around 2 x 10 19 Mx. 
4-3.2. Size of Elementary Flux Emergence 

We also measured the typical size of elementary struc- 
tures of emergence d unlt defined as the distance between 
two footpoints of individual fi-loops at their emergence 
period. Some small footpoints of emerged loops might be 
finally transported to the border of supergranules by the 
local convection. Others cause cancellation between the 
opposite polarities and disappear, which enlarges the dis- 
tance between two footpoints of individual O-loops, (i.e. 
dunit)- Thus the d U nit varies with time. For the accu- 
rate descriptions, we adopted the d U nit at the epoch when 
the fi-loops were observed as Call H filaments (~10 min- 
utes after the start of the emergence; Otsuji et al. 2007). 
Figure 5d shows the scatter plot of d max and d uni t, which 
suggest that d unit takes the values in the range of 2-6 
Mm irrespective of d max . Thus elementary and unitary 
O-loops in any emerging flux region has a common size of 
around 4 Mm, which is consistent with the most unstable 
wavelength (2-4 Mm) of Parker instability and preceding 
observation/simulation studies (Magara 2001; Pariat et al. 
2004; Isobe et al. 2007). 
4-3.3. Magnetic Flux Evolution 

The relation between the maximum flux $ max and the 
flux growth rate (d$/dt) is shown in Figure 5e. In the 
scatter plot, the data points distribute along the relation 
of 

9.6xl0 7 $ m 5 ax , (3) 



dt 



T 



where T is the emergence duration, (d$/dt) is in Mx hr _1 
and $ m ax in Mx. From equation (3), we can derive the 
emergence duration T in the unit of hour as a function of 
$max, which is 

T=1.03xKr 8 $ max . (4) 

Equation (4) indicates that an emergence event with 
large maximum flux shows relatively rapid magnetic flux 
growth. According to equation (4), when a flux tube of 
$o emerges with To, the flux tube with 2$o emerges with 
Ti = 2°- 43 T w 1.4T . T docs not depend linearly but non- 
linear ly on 'I'max- A tube with more magnetic flux emerges 
with less time than in the case of linear dependency. While 
the equation (4) is consistent with previous observations 
such as Zwaan (1987) and Hagennar (2001), the empirical 
relation (4) was first derived with wide range of magnetic 
parameters by Hinode high-resolution samples. 
4-3.4- Relation between Footpoints Separating Speed and 
Maximum Flux 
Figure 5f presents the relation between 'I'max and (i>), 
which indicates that the larger size EFRs show the sep- 



arating speed less than 1 km s _1 , while the small scale 
EARs footpoints separate with various speed up to 4 km 
s _1 . The mean separating speed (v) can be written as 

( W H^7-- (5) 



(G) 



From equations (2) and (4), equation (5) reduces to 
(v) = 2.1 x 10 3 $~ 016 , 

\ u / A " max ' 

where (v) is in km s _1 and $ ma x in Mx. Equation (6) 
indicates that the footpoints of emerged flux tube with 
less magnetic flux separate each other with larger speed. 
This relation is plotted with solid line in Figure 5e, which 
is consistent with the observed values. 



5. Discussion 

5.1. Size and Flux Dependence of SZ Formation 

The SZ features are associated with the magnetic flux 
emergence of more than around 2 x 10 19 Mx. Magara 
and Longcope (2003) performed three-dimensional MHD 
simulation of emerging magnetic flux and suggested that 
emerging field lines take the evolutionary path of a simple 
expansion if they emerge with a large aspect ratio (the ra- 
tio of their height to their footpoint distance); otherwise, 
field lines are inhibited from expanding and they show an 
undulating behavior (i.e. SZ formation). Our results pro- 
vide a new criterion of the total magnetic flux regarding 
the formation of SZs. 

5.2. Flux Dependence of the Spatial Size of Flux 
Emergences 

Another notable result is that the flux emergence phe- 
nomena with/without SZ follow the relation between the 
total flux and the maximum spatial size as described in 
equation (2). This relation is derived from wider range of 
magnetic parameters than the previous studies (Bruzek 
1967; Zwaan 1987; Harvey and Martin 1973; Harvey et al. 
1975). There authors gave the result for only narrow range 
of magnetic parameters. We present a comprehensive re- 
sult on this relation compared to the previous works. 

Let us try to derive the power-law relation between the 
maximum flux $ ma x and the maximum separation dis- 
tance d max , 



d cx <J> Ql 

"max ^ ^ ma; 



(7) 

from the viewpoint of simple dimensional analysis. Figure 
6 shows the schematic image of flux emergence from the 
convection layer. First, e? max is estimated as follows. 
Initially horizontal flux tube in the convection layer will 
rise with the typical length A = 10 ~ 20ifp where Hp is the 
local pressure scale height. The maximum separation dis- 
tance <i max between two main spots depends on the most 
unstable wavelength of Parker instability A at the initial 
depth of the flux tube, 

d max cxXcxH P . (8) 

Thus <i max is proportional to Hp at the depth where ini- 
tially the flux tube locates. 
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Large EFR 



Large flux tube 



Small flux tube 



Small EAR 




Fig. 6. Schematic image of emerging flux tubes, d is the 
distance between the main spots. I is the mixing length. Hp 
is the pressure scale height at the initial depth. B is the 
field strength in the flux tube, p is the mass density around 
the flux tube. A is the most unstable wavelength of Parker 
instability. 

Next the total flux 3> max is estimated as follows. From 
equipartition arguments, the magnetic and kinetic energy 
in the solar convection layer will balance to each other, 



8tt 



1 2 

2 P^conv 7 



where B, p and v conv are the field strength inside the 
flux tube, the mass density around the flux tube and 
the mean convection velocity. From mixing length the- 
ory (Stix 1989), v is given as 

fconv OC \/#p- ( 10 ) 

From equations (9) and (10), 

Bocv^. (11) 

Now we assume that the solar convection layer can be 
approximated by an adiabatically stratified atmosphere 
(Foukal 2004), 

Toc^ 1 . (12) 

T and 7 are temperature and adiabatic index 7 = cp/cy, 
where cp and cy are the specific heats at constant pressure 
and volume, respectively. The local scale height Hp is 
proportional to the temperature T, thus 



p oc T-1- 1 oc Hp 1 . 
From equation (11) and (13), 



(13) 



2(7-1) 



(14) 

Let us think about the flux tube width w. If the w is much 
larger or smaller than the local mixing length I oc Hp , the 



Photosphere 




Convection layer 



Fig. 7. Schematic image of flux tube underneath solar pho- 
tosphere. The gray areas represent the cross-section of flux 
tubes with the field strength of B. w and h are horizontal 
and vertical width of the flux tube, respectively, v is the rise 
velocity of the flux tubes at the photosphere. Small circles 
in the flux tube represent unitary and elementary flux tube 
found in Section 4.3.2. 

flux tube will be disintegrated by the convection flows or 
accumulated at the convection boundary. Thus the flux 
tube width is expected to be comparable to the mixing 
length, 



w ~ I oc Hp . 
Thus the total flux $ max can be estimated as 



w z B oc H 



2(7-1) 



(9) From equation (8) and (16), 



(15) 



(16) 



(17) 



Thus the power-law index of the relation between the 
maximum flux $ max and the maximum separation dis- 
tance (i m ax is derived to be ol\ = 2 g 7 ~^ . For example, 
a.\ is 0.30 with the adiabatic index 7 = 5/3 (ideal gas 
case realised in deep convective layers). With 7^4/3 at 
near the solar surface where the ionization status is chang- 
ing rapidly (Bhatnagar and Livingston 2005), a.\ is 0.25. 
These calculated values arc comparable to the observed 
value ai = 0.27. From the argument above, we get the 
view of the emergence depicted as in Figure 6. Magnetic 
tubes of large flux are anchored in deep layers and appear 
with large separation between two main spots on the solar 
photosphere. 

5.3. Flux Growth Rate 

We derived the relation between the maximum flux 
■JVax and the flux growth rate (d$/d£) to be equation 
(3). 

By a simple model of an emerging tube with uniform 
magnetic flux density, let us try to derive the power-law 
relation between the maximum flux and the flux growth 
rate, 



(d*/dt)oc*S 



(18) 
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Figure 7 shows the schematic image of the flux tubes 
just beneath the photosphere. Around the photosphere, 
the plasma (3 is almost 1. Thus the magnetic pressure 
B 2 /8ir ~ P = const, where P is the gas pressure at the pho- 
tosphere. So we can think that flux densities are nearly 
constant irrespective of the spatial size or total flux of the 
magnetic tube, 

B = const. (19) 

The rise velocity v of the flux tube is estimated as fol- 
lows. When the apex of flux tube reaches to the under- 
neath the solar surface, the rise motion is suppressed and 
the tube top becomes flattened (Magara 2001). Our obser- 
vation (Section 4.3.2) showed that the emergence occurs 
in unitary form irrespective of total magnetic flux. In 
these situations the rise velocity v from the photosphere 
will not depend upon the total magnetic flux of the tube, 

v = const. (20) 

Assuming that the aspect ratio of flux tube width h/w is 
constant just beneath the photosphere, 

h . . 

— = const. (21) 

w 

where w and h are horizontal and vertical width of the 
flux tube, respectively. 

The flux growth rate and the total flux are described as 

(d$/df) =wvB(xw, (22) 

^max = whB oc w h (X w 2 , (23) 

respectively. From equation (22) and (23), the relation 
between $ max and (d$/dt) are written as 

(d$/dt) oc $| ax . (24) 

Thus the power-law index is 

a 2 = 0.5, (25) 

which is consistent with the observed value 0.57. 

In this discussion we have not considered the factors 
such as magnetic field stratification inside the flux tube, 
realistic aspect ratio of tube and so on. If we include 
these factors then we can get more realistic interpretation 
of equation (3). 

6. Summary 

We investigated the morphological, dynamical and mag- 
netical characteristics of various flux emergence phenom- 
ena using high-resolution Hinode SOT data. To estimate 
the magnetic field density of SOT data we used SOHO 
MDI magnctogram data for the calibration. From 101 
samples of flux emergence events, we derived the total 
flux, flux growth rate, maximum separation and mean 
separation speed. The SZ features are associated with the 
magnetic flux emergence of more than around 2 x 10 19 Mx. 
The magnetic flux growth rate, emergence duration and 
mean separation speed were found to follow the power- 
law functions of the total magnetic flux with the indices 
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of 0.57, 0.43 and -0.16, respectively. The typical size of 
elementary emergence structures is around 4 Mm, which 
is consistent with the most unstable wavelength (2-4 Mm) 
of Parker instability. The mean separating speed (v) de- 
creases with larger magnetic flux. 

We got a physical view of solar flux emergence that 
emerging magnetic fields float and evolve balancing to 
the surrounding turbulent atmosphere from the discus- 
sion on the observed power-law relations. These observa- 
tional results should be verified by future numerical stud- 
ies. Possible influence of twisting or pre-existing magnetic 
field could be studied with the data of horizontal magnetic 
field in the emerging site and will be reported in the near 
future. 
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# 


Date 


^max 


(v) 


$max 


(d$/dt) 


SZ 


# 


Date 


^max 


(v) 




(d$/di) 


SZ 


1 


20061126 


15.2 


0.97 


1.2el9 


2.6el8 


N 


51 


20080807 


7.4 


1.87 






N 


2 


20061201 


71.0 


0.62 


3.0e21 


9.5el9 


Y 


52 


20080809 


3.9 


3.09 





— 


N 


3 


20061209 


16.6 


0.36 


9.2el9 


7.1el8 


Y 


53 


20081014 


12.2 


1.07 





— 


Y 


4 


20061226 


24.9 


0.74 


— 


— 


Y 


54 


20081018 


7.5 


1.20 


3.9el8 


2.3el8 


N 


5 


20070104 


17.4 


0.58 


2.0c20 


2.4el9 


Y 


55 


20081019 


8.4 


1.09 


6.5cl8 


3.0el8 


N 


6 


20070117 


17.3 


0.26 


— 


— 


Y 


56 


20090118 


14.3 


0.73 


1.7e20 


3.2cl9 


Y 


7 


20070204 


10.7 


0.33 


1.2e20 


1.3el9 


Y 


57 


20090226 


17.6 


0.56 


2.6el9 


3.0el8 


N 


8 


20070205 


19.2 


1.21 


1.5e20 


3.3el9 


Y 


58 


20090314 


10.9 


1.45 


— 


— 


N 


9 


20070213 


7.2 


1.47 


5.0el8 


3.7el8 


N 


59 


20090402 


8.8 


0.86 


3.6el9 


1.3el9 


N 


10 


20070219 


8.5 


2.72 


3.6el8 


4.1el8 


N 


60 


20090426 


18.5 


1.44 


1.2el9 


3.4el8 


Y 


11 


20070308 


16.8 


1.38 


5.7el9 


1.7el9 


Y 


61 


20090601 


32.6 


2.48 


2.7e20 


7.3cl9 


Y 


12 


20070328 


7.8 


0.29 


— 


— 


N 


62 


20090603 


30.7 


0.54 


3.1e20 


2.0cl9 


Y 


13 


20070331 


11.6 


1.76 


6.6el8 


3.6el8 


N 


63 


20090623 


7.4 


2.52 


l.lel9 


1.3cl9 


N 


14 


20070414 


10.7 


1.70 


2.1cl9 


1.2el9 


Y 


64 


20090704 


63.4 


0.35 


— 


— 


Y 


15 


20070419 


16.7 


1.09 


4.5cl9 


1.0el9 


Y 


65 


20090707 


51.6 


0.54 


2.3e21 


8.8cl9 


Y 


16 


20070601 


15.3 


0.50 


— 


— 


Y 


66 


20090821 


5.9 


1.01 


4.6cl8 


2.8cl8 


N 


17 


20070603 


14.2 


0.73 


— 


— 


Y 


67 


20090928 


10.2 


1.18 


2.8el9 


1.2el9 


N 


18 


20070702 


8.2 


0.93 


1.8el9 


7.2el8 


N 


68 


20091009 


9.8 


1.47 


2.3el8 


1.2el8 


N 


19 


20070807 


32.0 


0.46 


4.4e20 


2.3el9 


Y 


69 


20091015 


9.4 


2.09 


1.2cl9 


9.6cl8 


N 


20 


20070825 


12.9 


1.90 


4.2el9 


2.2el9 


Y 


70 


20091019 


10.0 


1.29 


4.5cl8 


7.8cl8 


N 


21 


20070918 


9.6 


2.49 


9.7cl8 


1.8el9 


N 


71 


20091023 


34.9 


0.44 


— 


— 


Y 


22 


20070929 


20.6 


3.99 


1.3el9 


9.4el8 


Y 


72 


20091025 


9.7 


1.00 


— 


— 


N 


23 


20071001 


22.5 


1.48 


l.le20 


2.7el9 


Y 


73 


20091026 


49.8 


1.91 


1.7e21 


2.4e20 


Y 


24 


20071008 


11.2 


1.70 


2.5el9 


1.4el9 


Y 


74 


20091207 


11.2 


1.06 


6.0cl8 


2.1el8 


N 


25 


20071009 


16.3 


1.42 


2.0el9 


6.4el8 


Y 


75 


20091215 


80.5 


2.00 


9.0e20 


8.0el9 


Y 


2G 


20071015 


6.4 


0.59 


— 


— 


N 


76 


20091226 


19.7 


1.03 


1.4e20 


2.7el9 


Y 


27 


20071020 


8.8 


0.63 


— 


— 


N 


77 


20091229 


32.6 


1.03 


5.5e20 


6.3el9 


Y 


28 


20071026 


8.7 


1.36 


1.4el9 


7.7cl8 


N 


78 


20091230 


67.3 


0.57 


3.0e21 


9.1cl9 


Y 


29 


20071108 


4.4 


3.69 


2.1el8 


6.4el8 


N 


79 


20100108 


28.4 


0.67 


6.0e20 


5.1cl9 


Y 


30 


20071111 


19.1 


3.35 


1.6el9 


1.0el9 


Y 


80 


20100109 


45.0 


0.69 


2.6e21 


1.4e20 


Y 


31 


20071117 


7.4 


1.41 


1.2cl9 


8.1cl8 


N 


81 


20100110 


95.6 


0.57 


6.7e21 


1.4e20 


Y 


32 


20071123 


24.8 


0.63 


— 


— 


Y 


82 


20100124 


15.7 


0.94 


1.3e20 


2.8el9 


Y 


33 


20071127 


12.5 


2.51 


4.3el8 


3.1el8 


N 


83 


20100217 


25.3 


0.71 


1.3e20 


1.3el9 


Y 


34 


20071129 


13.2 


1.11 


— 


— 


Y 


84 


20100221 


10.7 


1.75 


7.9el9 


4.7el9 


Y 


35 


20071207 


13.2 


0.33 


— 


— 


N 


85 


20100222 


32.5 


2.20 


3.3el9 


8.1el8 


Y 


36 


20071208 


35.6 


0.65 


1.2e21 


8.1cl9 


Y 


86 


20100313 


11.6 


1.22 


2.4el8 


9.0el7 


N 


37 


20071209 


17.2 


0.96 


9.6el9 


1.9el9 


Y 


87 


20100322 


4.5 


1.17 


1.8cl8 


1.6cl8 


N 


38 


20071210 


24.9 


1.56 


1.6e20 


3.6el9 


Y 


88 


20100325 


11.3 


1.46 


3.8cl9 


1.8cl9 


Y 


39 


20071211 


27.1 


0.59 


9.5e20 


7.5el9 


Y 


89 


20100414 


10.3 


1.28 


3.1cl9 


1.4el9 


N 


40 


20071212 


15.4 


0.75 


— 


— 


Y 


90 


20100612 


46.0 


0.19 


l.le20 


7.2el9 


Y 


41 


20071213 


24.7 


0.66 


6.3e20 


6.0el9 


Y 


91 


20100613 


28.1 


0.65 


— 


— 


Y 


42 


20080105 


11.9 


0.45 


— 


— 


N 


92 


20100618 


13.8 


0.51 


5.8el9 


7.7el8 


Y 


43 


20080106 


13.9 


1.00 






N 


93 


20100619 


38.0 


0.55 


l.le21 


5.9el9 


Y 


44 


20080113 


18.1 


1.09 






Y 


94 


20100723 


18.5 


1.03 


2.4e20 


4.8cl9 


Y 


45 


20080124 


18.2 


1.30 


3.0cl9 


7.8el8 


Y 


95 


20100727 


13.6 


0.80 


6.0cl9 


1.3cl9 


Y 


4G 


20080419 


17.1 


1.41 


7.7el9 


2.3el9 


Y 


96 


20100728 


48.4 


0.56 


1.8c21 


7.6cl9 


Y 


47 


20080610 


21.8 


0.84 


8.3el8 


l.lel8 


N 


97 


20100805 


28.4 


0.42 


5.2e20 


2.8el9 


Y 


48 


20080615 


16.2 


0.33 






Y 


98 


20100809 


10.6 


0.27 


1.7e20 


1.5el9 


Y 


49 


20080620 


4.6 


2.08 






N 


99 


20100812 


11.6 


0.29 


7.1el7 


1.0el8 


N 


■50 


20080630 


3.6 


1.56 


2.6el7 


4.1el7 


N 


100 


20100814 


43.9 


1.03 


4.2e20 


3.5el9 


Y 



101 20100822 12.4 2.75 7.9el8 6.3cl8 N 

Column(l): identification number of the flux emergence. Column(2): date of the flux emergence. Column(3): maximum 
distance between two main spots (Mm). Column(4): mean separating speed of the main spots (km s" 1 ). Column(5): maximum 
flux increment (Mx). Column(6): flux growth rate (Mx hr -1 ). Column(7): existence of stagnation zone (Yes/No). 
— The missing value ( — ) is due to the event with only Call H observation and without magnetogram. 



12 



K. Otsuji et al. 



Ichimoto, K., et al. 2004, SPIE, 5487, 1142 

Isobe, H., Tripathi, D., & Archontis, V. 2007, ApJ, 657L, 53 

Kitai, R. 1983 Sol. Phys., 87, 135 

Kosugi, T., et al. 2007, Sol. Phys., 243, 3 

Kurokawa, H., Kawaguchi, I., Funakoshi, Y., & Nakai, Y. 1982, 

Sol. Phys., 79, 77 
Lamb, D. A., DeForest, C. E., Hagenaar, H. J., Parnell, C. E., 

& Welsch, B. T. 2010, ApJ, 720, 1405 
Loughhead, R. E., & Bray, R. 1961, Australian Journal of 

Physics, 14, 347 
Magara, T. 2001, ApJ, 549, 608 

Magara, T. & Longcope, D. W. 2001, ApJ, 559, L55 
Magara, T. & Longcope, D. W. 2003, ApJ, 586, 630 
Martin, Sara F. 1988, Sol. Phys., 117, 243 
Matsumoto, R., & Shibata, K. 1992, PAS J, 44, 167 
Matsumoto, R., Tajima, T., Shibata, K., & Kaisig, M. 1993, 
ApJ, 414, 357 

Matsumoto, T., Kitai, R., Shibata, K., Otsuji, K., Naruse, T., 

Shiota, D., & Takasaki, H. 2008a, PASJ, 60, 95 
Matsumoto, T., et al. 2008b, PASJ, 60, 577 
Matsumoto, T., & Kitai, R. 2010, ApJ, 716, L19 
Murray, M. J., et al. 2006, A&A, 460, 909 
November, L. J., & Simon, G. W. 1988, ApJ, 333, 427 
Nozawa, S., Shibata, K., Matsumoto, R., Sterling, A. C, 
Tajima, T., Uchida, Y., Ferrari, A., & Rosner, R. 1992, 
ApJS, 78, 267 
Nozawa, S. 2005, PASJ, 57, 995 
Otsuji, K., et al. 2007, PASJ, 59, 649 

Otsuji, K., Kitai, R., Matsumoto, T., Ichimoto, K., Ueno, S., 

Nagata, S., Isobe, H., & Shibata, K. 2010, PASJ, 62, 893 
Pariat, E., Aulanier, G., Schmieder, B., Georgoulis, M. K., 

Rust, D. M., & Bernasconi, P. N. 2004, ApJ, 614, 1099 
Parker, E. N. 1966, ApJ, 145, 811 
Scherrer, P. H., et al. 1995, Sol. Phys., 162, 129 
Shibata, K., Tajima, T., Steinolfson, R.S., & Matsumoto, R. 

1989, ApJ, 345, 584 
Shimizu, T., et al. 2007, Sol. Phys., 249, 221 
Stix, M. 1989, Springer- Verlag, The Sun, 198 
Strous, L. H., Scharmer, G., Tarbell, T. D., Title, A. M., & 

Zwaan, C. 1996, A&A, 306, 947 
Strous, L. H., & Zwaan, C. 1999, ApJ, 527, 435 
Suematsu, Y., et al. 2007, Sol. Phys., 249, 197 
Titov, V. S., Priest, E. R., & Demoulin, P. 1993, A&A, 276, 

564 

Tsuneta, S., et al. 2007, Sol. Phys., 249, 167 

Watanabe, H., Kitai, R., Okamoto, K., Nishida, K., 

Kiyohara, J., Ueno, S., Hagino, M., Ishii, T. T., & 

Shibata, K. 2008, ApJ, 684, 736 
Zirin, H. 1972, Sol. Phys., 22, 34 

Zwaan, C, Brants, J. J., & Cram, L. E. 1985, Sol. Phys., 95, 
3 

Zwaan, C. 1987, Annual review of astronomy and astrophysics, 
25, 83 



